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Abstract. When a large group of pedestrians moves around a corner, most pedes- 
trians do not follow the shortest path, which is to stay as close as possible to the inner 
wall, but try to minimize the travel time. For this they accept to move on a longer 
path with some distance to the corner, to avoid large densities and by this succeed in 
maintaining a comparatively high speed. In many models of pedestrian dynamics the 
basic rule of motion is often either "move as far as possible toward the destination" or 
- reformulated - "of all coordinates accessible in this time step move to the one with 
the smallest distance to the destination. Atop of this rule modifications are placed to 
make the motion more realistic. These modifications usually focus on local behavior 
and neglect long-ranged effects. Compared to real pedestrians this leads to agents in a 
simulation valuing the shortest path a lot better than the quickest. So, in a situation 
as the movement of a large crowd around a corner, one needs an additional element 
in a model of pedestrian dynamics that makes the agents deviate from the rule of the 
shortest path. In this work it is shown, how this can be achieved by using a flood fill 
dynamic potential field method, where during the filling process the value of a field 
cell is not increased by 1, but by a larger value, if it is occupied by an agent. This 
idea may be an obvious one, however, the tricky part - and therefore in a strict sense 
the contribution of this work - is a) to minimize unrealistic artifacts, as naive flood fill 
metrics deviate considerably from the Euclidean metric and in this respect yield large 
errors, b) do this with limited computational effort, and c) keep agents' movement at 
very low densities unaltered. 
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1. Introduction 

1.1. Pedestrians in the Light of Light 

Everyone knows from own experience that even when selecting the mode of travel but 
also as a participant in any mode of traffic minimizing the travel time is an important 
criterion when planning a trip. There can be other criteria against which the shortest 
travel time is weighed against for example monetary costs or sometimes pedestrians 
value staying away from large streets high. However, minimizing the travel time is 
almost never irrelevant. If the travel time is not minimized for a full path one can often 
insert intermediate destinations on a path, to make the rule of minimal travel times 
apply for each of the partial paths. 

As a first step of modeling let's neglect all other criteria and assume for pedestrians: 

Pedestrians try to move from their current position to their destination on the 
path of shortest travel time. 

This bears a remarkable similarity to the formulation of Fermat's Principle for the 
path of light: 

On the way from a given point A to a given point B a ray of light takes the 
path of least izme.jl| 

Why then is the motion of pedestrians and light so different? It's not only the kind 
of interaction, but a more fundamental difference: Fermat's Principle is only seemingly 
teleological, while in fact it is the result of laws of nature following causality. Classically 
this follows from Huygen's Principle ^ and quantum-mechanically from Feynman's 
Path Integral Formalism I.e. for certain calculations the teleological and the causal 
formulation are equivalent. Pedestrians on the contrary act teleologically as they have 
- at least in introspective - a free will and make plans. Additionally causal restrictions 
apply to their motion. If one restricts the motion of participants in an experiment to 
a one-dimensional motion, the focus is set to these causal restrictions. In such a case 
there is no more teleological element present. An experiment purely investigating this 
is for example reported about in 0, 0]- In other cases experiments and observations 
often do not exclusively investigate the causal restrictions but still emphasize it a lot, 
as there is usually a very limited set of alternatives available for any kind of planning 



As a result of this double influence - causal and teleological - and because of limited 
access to information of the path ahead - pedestrians only try to minimize their time 



process |6l4l3 . 
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and usually only partially succeed with that. On the contrary for light there are only few 
clearly defined and well understood exceptions from the path of minimal time: in very 



intense (many photons) ray of light with very high energy^ 



non-negligible probability that light scatters at light [ij] and the ray is "fanned out" 



frequency) it can happen with 



and - just for completeness - gravitational lenses can produce images with light moving 



on a path of (locally) maximal time [15| . The case of gravitational lenses is not of much 
interest here. But light-light scattering is of some interest as with it the applicability 
of Fermat's Principle vanishes for light just for a similar reason why Fermat's Principle 
cannot be used to do calculations of pedestrian dynamics. Within a normal course of 
life beam of light photons typically do not interact as the photon-photon scattering 
probability (its cross section) is by far too small. 

As long as the photons within a ray of light do not scatter, Fermat's Principle 
expressed in a formula says that the travel time 

r = ^ (1) 

J PathP:x=A '^i^) 

is minimal for the path P "chosen" by the ray of light. v{x) = c/n{x) is the speed 
of light as the case may be within some material with some refractivity n{x). From 
knowing the field of speeds \v{x)\ one could calculate the field T{x) of minimal travel 
times to the destination B. For each position the direction of the ray of light then would 
be equal to the one of the gradient VT(x), as Fermat's principle compulsorily also holds 
for each partial path. 

Classical optics is simpler than high energy quantum optics or the movement of 
pedestrians as v{x) only depends on the material and its refractivity at position x and 
not the light itself or to be more precise the spatial-temporal distribution of photons. 
For pedestrians v{x, t) heavily depends on the local distribution of other pedestrians at 
a certain time. Put mathematically: the velocity field v{x, t) - or even more general: 
the field of velocity distributions - is not an input but a result of the calculation of 
pedestrian dynamics and those space-time regions are of specific interest, where the 
speed is significantly lower than the desired speeds: "jams". Obviously one is faced 
with a highly non-linear phenomenon. 

Having realized that the most important input factor is not available in pedestrian 
dynamics, the first and radical step is to assume that the agents do not influence each 
other at all. This is equivalent to allowing two agents to take the same place at the 
same time, which of course is not realistic and would result in all agents walking the 
same spatially shortest path regardless of the density on that path. So, the minimal 
mutual influence of agents is a strong, short-ranged, repulsive influence, often simply a 
hard-core exclusion, respectively a hard-core repulsion. 

I The reason for this being that charge conservation does not allow a tree-level process. So the 
scattering in lowest order involves two electron propagators which in the end lead to the electron 
mass to the eighth in the denominator. Furthermore the four vertices bring in an a**, which decreases 
the scattering probability further. Only if energy and density of the photons is large enough these 
factors can be balanced to non-negligible overall scattering probabilities. 
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As apart from this exclusion principle at this stage no additional interaction between 
agents exists, the basic direction of motion then is determined from the gradient of the 
field of spatial distances to the destination: VS'(x). The element of a model that 
incorporates the influence of VS'(x) mathematically often has the same structure as the 
other elements of the model. Its special status is then reflected in the typically most 
important role the parameters configuring the strength of the influence of V5'(x) have. 

Then as a third step the actual modeling process sets in, where all the available 
models differ from each other and deviations from the basic direction of motion are 
calculated. Normally at least one of these additional elements is meant to compensate 
for the agents following VS'(x) too strictly, i.e. walking too close to the spatially shortest 
path and therefore producing unrealistic jams. At this point the modeler is faced with 
a problem: In many situations this element of the model would have to incorporate 
quite long-ranged interactions to make the movement of agents realistic. But long- 
ranged interactions often are computationally costly. In a majority of cases modelers 
decide in favor of quick computation and against long-ranged interactions, accepting the 



formation of unrealistic jams in certain situations 16|, as the agents esteem the spatially 



shortest path too much compared to the quickest path. For a notable exception, where 



a minimized travel time is considered in an integrated way, see |17l . 

This work investigates the question, if the field (aka "potential") S{x) can be altered 
with some method to some field S{x, t) (here a time dependence must be introduced, 
because S{x,t) necessarily depends on the distribution of agents), such that the basic 
direction of motion in general is more similar to the if it had been calculated 

from the unknown T{x,t). 

It is assumed that the pedestrians to be simulated have fairly good knowledge of 
their environment. Contrary to an ansatz like in [l9], where agents' poor knowledge of 
the environment is intended to be handled in the operational model, this aspect is spared 
here, as it is believed that this should be handled on the tactical or strategical level of 
scenario modeling: an agent has perfect knowledge of the spatial structure relevant 
for the path between his current origin-destination relation. If this was not the case, 
the scenario modeler would have to split up the origin-destination relation into smaller 
pieces until the condition holds. Incomplete information then would be reflected in the 
routing decision at an origin with multiple potential destinations attached and not all 
over the path between origin and destination. For a description of the different levels of 



modeling see for example 20 



1.2. The Use of Potentials in Pedestrian Dynamics and Robotics 



Potential fields have long been used in robotics 2l|, |22] (let's call them U{x) here) and 



the simulation of pedestrians |23l-l29l| (let's call them S{x) here) to control robots and 
agents to evade obstacles and move toward a destination. The difference is that in 
robotics it is usually assumed that an autonomous robot knows about its destination 
coordinate but has no knowledge of the position of obstacles except for those which it 
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"sees" . For pedestrians on the contrary it is typically assumed that they have at least 
some knowledge on the whole path, the positions of static obstacles and the detours 
they have to walk compared to linear distance, even if they actually only see a small 
fraction of the whole path. This has consequences for the calculation and use of the 
potential. In robotics the artificial potential UaX position x from a destination at Xd 



was originally [2l| calculated as 



U artificial (•^) 



(2) 
(3) 



In such a potential local minima can occur and a robot has to be equipped with the 
ability to realize it is in a minimum and how to get out of it. Much different from this is 
the basic assumption in many models of pedestrian motion that pedestrians have a good 
global knowledge of the exact Euclidean distances and know about the shortest path 
between their current position and their destination under consideration of all obstacles. 
This assumption can be unrealistic for very complex geometries like huge mazes, but 
for many situations it comes close to reality. In this contribution "potential" is used in 
the sense of pedestrian dynamics. 

The double usage of the word "potential" has repeatedly led to confusion. For 
example in 19|] a number of models have been attributed to use robotics-like potentials, 
although they do not. In contrast the "sink propagation value" method proposed there 



is very similar to the static floor field used in 30| rendered more precisely in 26|, and 



31 



In publications of the Cologne- Tokyo school of pedestrian dynamics modeling 
23l . |26[ | the potential is called "static floor field" . Similar to 19| the values of the static 



floor field decrease with increasing distance from the destination, linearly reflecting the 
attractiveness of a cell. In [32i] a computer science perspective has been adopted: the 
distance to the destination is a relevant input value, so a look up table of distances 
is written to memory. The difference to the static floor field is trivial, but to not get 
messed up and miss a sign, the latter field has been called "distance potential field" . 

There are many methods to calculate a potential S{x) that consider obstacles, when 
calculating distances. Flood fill methods, Diikstra's algorithm on a visibility graph 
26l . l28l . I33I . l34| . the fast marching method [35| and ray casting |32l |. These methods 
have been compared in [32] or used and applicated in [36| and it has been concluded 
that in large real world applications only flood fill methods are quick enough to be 
calculated each time step or each N time steps, with N chosen such that the calculation 
is done for example each simulation second. Flood fill methods, however, yield large 
errors compared to true Euclidean distances [32 . 



1.3. Floor Field Models of Pedestrian Dynamics 



In floor field models 2314281. l30l . l3ll . l37l . l38l ] agents move on a discrete lattice from 



cell to cell. Before an agent moves, a probability is assigned to each cell he could 
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reach and then a destination cell is selected randomly according to these probabilities. 
Most influences on an agent's motion are modeled as partial probabilities following an 
exponential function 

= (,kx-func{X)^ 

where kx is a coupling constant, denoting the strength of the influence. The partial 
probabilities are multiplied to form the full probability that a cell is chosen as destination 

P=N\{px. (5) 

X 

N normalizes the probability such that the sum over all accessible cells equals 1. For 
this work it is sufficient to assume that all previously introduced influences, except for 
the influence of the static floor field (resp. distance potential field) - containing in each 
cell the information on the distance of that cell (x, y) to the destination - are switched 
off: fcx = 0, except for kg, 

p,y = Ne-^^'^^y (6) 



2. Altering the Distance Potential Field — A Quick Method as Perturbation 
of an Exact One 

As has been stated, there are methods that are sufficiently precise in the calculation 
of distances, such that no unrealistic looking artifacts in the motion of the agents can 
be observed in a simulation. But these methods need painfully large times for their 
calculation in large scenarios, if they are executed each simulation second. And there 
are methods, which are in principle quick enough to calculate the potentials each time 
step, but these methods often make the agents move on strange trajectories, as they are 
not very precise with regard to Euclidean metric. 

The idea is, to calculate an exact and precise {S = S{x)) potential and a potential 
using a quick method {Sdyn = Sdyn{x,t = 0)) potential in advance, i.e. before any agent 
has entered the simulation. And then calculate each time step - dynamically - with the 
quick method a potential under consideration of the agents: = Sdyn{x,t)). The 

probability to select a cell as destination is then calculated as 

p = Are-'=s(s+Sd,„-5S,n) (7) 

For sake of clarity the x and y indices of equation have been dropped. With this 
method the imprecise but quick method becomes a perturbation of the exact method. 
If there is only one agent in the scenario, he is only minorly affected, as all perturbation 
of the field is "behind" him, i.e. further away from the destination than he himself is. 
As only differences between S^^^y^ and S^y^ influence the motion, one can hope that this 
difference also cancels out a large share of the errors of the quick method. Introducing 
a new coupling constant for the new effect, one can generalize equation ([71) to 

p = Are-'=s^-'=s(-5d,„-5S,„) = A^e-''s^e"''s('5d.n-5S,„) = NpsPs (8) 
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In this way an entirely new influence has been created, fully independent from the 
influence of the static floor fleld. It cannot be called dynamic floor fleld, as this is the 
name for the virtual trace 23|, 1261], so it is referred to as "dynamic distance potential 
fleld" here. 



3. The Quick Method: Flood Filling 

It is not only that the exact methods take too long for their calculation, also only a 
subset of them would be able to consider the agents in one way or another meaningfully 
as a modiflcation. This is different with flood fill methods. In the process of flooding 
the scenery, one simply adds 1, if the cell is not occupied and a value Sadd > 1, if it 
is occupied by an agent. Basically there are two simple flood flU methods: the flrst 
acts in the von Neumann neighborhood. It is called the Manhattan metric and is the 
p = 1 vector norm. The other one acts on the Moore neighborhood. It is called the 
Chebyshev metric and is the p — oo vector norm. A combination of these is variant 1 



(VI) as introduced in [32], where the maximum error is limited to smaller values than 
in the two basic methods on cost of having to do two flood flUs and calculate a square 
root for each cell. 

The method proposed introduces two new parameters to the model: the coupling 
ksdyn, and the summand Sadd in case a cell is occupied. 



3.1. The Role of Parameter Sadd 

While the effect of ksdyn is obvious and just the same as the effect of all other coupling 
constants in a floor fleld model, the situation is a bit different for Sadd- A large Sadd 
does not in general increase the strength of the effect, this is only true for those cells on 
which agents are located. For all other cells Sadd determines the size a block of occupied 
cells needs to have that the flood partially flows through the block and not around. This 
does not only depend on Sadd but also on the shape of the block: there's only an effect, 
if from the potential's wave front perspective the block is wider than deep. If three 
agents are standing next to each other, all on the same value of the distance potential 
field, the value behind the central one is 3, if Sadd = 2 and 4 for any other value of Sadd 
(assumed Sadd £ A/" and Manhattan metric is used). If the solid block of occupied cells 
is deeper than wide, the flood always has to flow around the block. In other words: for 
large values of Sadd blocks of occupied cells act as if they were obstacles and the flood 
has to flow around them during the calculation of the potential. Only if such blocks 
would bar the flood's way to some part of the area or if Sadd is small and the blocks 
wider than deep the flood can also flow through these blocks. The former case can be 
interpreted as a block being inflnitely wide (from border to border). 

It's also blocks of agents ranging from border to border - not the idea that agents 
are in some way "softer" than obstacles - that are the reason, why agents can not be 
treated like obstacles, i.e. that Sadd is set to inflnity, respectively that the dynamic 
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distance potential field may not flood "over" agents at all. Otherwise the flooding 
algorithm could not reach parts of the area, if somewhere there is a chain of agents 
without a gap from one border to another. 

So there are two interpretations of Sadd'- if there is no line of agents ranging from 
one border to another Sadd is an estimation for the time to walk around the block, the 
bigger the block, the larger the time. If the block does range from one border to another 
Sadd is a measure for the speed reduction or time loss. The length of a jam ranging from 
border to border times Sadd — 1 then determines the length of an accepted detour. 



3.2. Examples 

Figures [H to [3] show the impact of a group of agents on S^y^ by showing S^^y^ 



_ cO 

dyn' 



It can be seen clearly, that with Manhattan metric the impact spreads infinitely in 
vertical (and horizontal) direction and with Chebyshev metric in diagonal direction. 
Only with VI the impact range is finite in any case. The seemingly strange pattern 
into diagonal direction with VI metric is a result of rounding effects: instead of all cells' 
value approaching to zero, with increasing distance less and less cells have a value 1. 








I 


•m 

\ 


1 











Figure 1. 5^,y„ - S'^y^ {sadd = 10; from left: Manhattan, Chebyshev, VI) with 17 
agents placed on two lines at the brightest white spots (exits are shown green, walls 
red). The difference of the time dependent and the empty potential is the larger, the 
brighter a spot is. 



4. Comparison to a Static Method 

The main problem, why this method is introduced is that without some method to 
prevent agents to crowd too densely at the inner side of the corner, the flow is reduced 
unrealistically at the corner, as only a tiny space at the corner is used by the agents. It 



has been shown [26|, |28| that this can be achieved partially by a simpler, static method, 
where the probability that a cell is chosen as destination is reduced, if its distance W 
to a wall is smaller than some distance Wmax'- 
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Figure 2. S^y^ - S^y^ [sadd = 10; from left: Manhattan, Chebyshev, VI) with 105 
agents placed at the bright white squares. The dilference of the time dependent and 
the empty potential is the larger, the brighter a spot is. 




However, this method either assumes that large amounts of agents will flow around a 
corner. Then Wmax has to be large to distribute the agents on a wide variety of radii 
from to Wmax around the corner. But then, if a small number of agents arrives, 
they will keep an unrealistically large distance to the corner. Or, if one expects small 
numbers, Wmax has to be small, and for large groups the effect will be too small. 

The stadium example of subsection 16.61 was calculated with the static method to 
make a comparison; see the discussion there. 



5. Transferability 



27 



28l . |37| as an example. 



The method was introduced above using the F.A.S.T. model 
but its application is by no means limited to this or any other floor fleld model. It 
can for example also be used in a continuous space model like the social force model 



13l . l36l . |39N42| . For a model that is continuous in space and not cell-based, one has to 



integrate the cell concept to the else continuous concept and therefore answer questions 
like "Should agents always be assigned to one cell, or can there be parts of agents on 
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cells?" , "What to do, if more than one agent occupies a cell?" , or "How large should the 
cells be chosen?" . Probably details of the answer will only affect details of the results. 
Finally, ^S\y^ = S\y^ — Sdyn can be used as a modification in the method, where the 
direction of the desired motion is determined. 

6. Discussion: Examples of the Influence of the Method 

6.1. The Influence on Free Speed 

The influence on free speed of this method is model specific. For the tests with the 
F.A.S.T. model the following settings were made: ks = 1.0, ksdyn = 10.0, kanyother = 0.0, 
and Sadd = 10. All simulations were repeated 100 times with identical settings. One 
single agent was placed 250 cells away from his destination. The room was 400 X 
400 cells wide, such that border effects could if at all only play an insignificant role. 
The following results for the number of rounds to reach the destination were obtained 
(format: average ± std (min, max)): 





none 


Manhattan 


Chebyshev 


VI 


1 


437.2 ± 21.6 


368.5 ± 14.0 


400.3 ± 19.9 


401.0 ± 21.2 




(388, 492) 


(326, 397) 


(346, 459) 


(351, 454) 


5 


61.32 ± 2.18 


60.87 ± 2.56 


60.97 ± 2.71 


61.15 ± 2.33 




(57, 69) 


(56, 67) 


(55, 69) 


(57, 67) 



So there is a some influence at Vmax = 1 and no recognizable influence at v^ax = 5 
and for Vmax = 1 the influence is stronger, when the Manhattan metric is used. This is 
due to the motion along the discretization axis. 

These results can be obtained analytically. If one sets the smallest value of the static 
distance potential field and the dynamic distance potential field to 0, it will always look 
like this, if the agent is meant to move downward: 
Exact static distance potential field: 



2 


2 


2 


1 


1 


1 












The Manhattan and Chebyshev dynamic distance potential fields are shown in table 

m 

This leads to exponents in the probabilities (with ks = 1.0 and ksdyn = 10.0) as 
shown in table [2] 
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2 


3 


2 


1 


10 


1 












2 


2 


2 


1 


10 


1 












Table 1. Manhattan and Chebyshev dynamic distance potential fields. 



-2 


-12 


-2 


-1 


-91 


-1 












-2 


-2 


-2 


-1 


-91 


-1 












Table 2. Exponents of probabilities for Manhattan and Chebyshev metric. 



and therefore the following expectation values for the time to reach the destination: 

r_e(250) = 250/ 3,0 + 3,-1 +3^-2 « 434.6 

, , ,1 -360 + 0(26-1 + 1 -6-91) + (-l)(2e-2 + e-i2) 

TManhattan{2hQ) = 250/ 3,0 + 26'! + 6-91 + 26-^ + 6-^2 ^ 366.9 

, , ,1 • 36° + 0(26-1 + 16-91) + (-1) • 36-2 

rc.e.,..e.(250) = 250/ 3^0 + 26-1 + 6-91+36-2 " ^^^-2 

If one does the same with Vmax = 5 one sees that all terms that bring in differences 
to the case without dynamic distance potential are small compared to those terms that 
are there in any case, which makes all changes vanish. Graphically the reason is that 
the neighborhood for Vmax = 5 is wider, but there is still at maximum one cell influenced 
by the dynamic distance potential. This in turn implies that if there are a few agents in 
the neighborhood of an agent, there can again be an accelerating effect at all maximum 
speeds. But all of this depends on the ratio ksdyn/ks- If it is small this effect will also be 
small and it must not be forgotten that for this investigation it was set to an absurdly 
large value of 10. 

6.2. A Crowd Moving Around a Corner 

This is the elementary situation for which the method was introduced, i.e. where the 
biggest problems occur, if the basic direction is calculated from S{x). As shown in figure 



Figure 4. Starting position of 1250 agents (red) to move along a corridor (black) 
around a corner to a destination (green). The walls are white. The inflow was restricted 
by a number of columns. 

m 1250 agents had to move around a corner. The maximum speed was set to 5 for all. 
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All simulation runs were repeated 100 times. Without dynamic distance potential field 
the process took between 308 and 319 rounds, at an average of 314.1 ± 2.4 rounds. The 
results for the various metrics and different parameters is shown in [3l 



ksdyn I Sadd 



Manhattan 



Chebyshev 



VI 



1.0 / 2 
3.0 / 2 
1.0 / 4 
3.0 / 4 
1.0 / 10 
10.0 / 10 



224.5 ± 2.2 (219, 231) 
191.4 ± 2.0 (186, 197) 
185.0 ± 1.7 (182, 190) 
178.7 ± 1.8 (174, 183) 
176.0 ± 2.1 (171, 181) 
190.3 ± 2.5 (185, 198) 



267.8 ± 1.8 (263, 267) 
242.1 ± 1.9 (237, 247) 
246.7 ± 1.8 (243, 252) 
235.7 ± 2.0 (231, 241) 

243.9 ± 1.7 (240, 249) 
240.5 ± 2.5 (234, 247) 



249.5 ± 1.9 (246, 255) 

214.1 ± 1.8 (211, 220) 

207.3 ± 2.0 (203, 215) 
181.7 ± 1.5 (178, 185) 

184.2 ± 1.7 (181, 190) 

196.4 ± 2.2 (190, 205) 



Table 3. Number of rounds until the last agent has walked around the corner and left 
the simulation 



Without dynamic potential after 90 seconds the situation looks like fig [5] 




Figure 5. Situation after 90 seconds without the dynamic distance potential field 
method. 



Figure [6] shows that the influence on the shape of the crowd is strongest if the 
Manhattan metric is used and weakest with the Chebyshev metric. This is not 
astonishing, as in this case the possibility to flood diagonally over just common corners 
and not only over common edges quickly can erase the effect of extra points in the 
potential due to agents. Furthermore one can see that large ksdyn tend to exhibit 
artifacts in form of unrealistic density waves and this in the end leads to a re-increase 
of the evacuation time, when ksdyn is too large. 

To check the direction dependence of the method, the geometry was rotated by 
45 degree. A direction dependency can come into play by the anisotropy of the speed 
neighborhood |43| and the anisotropy of the dynamic distance potential field method. 
Without dynamic distance potential field the process took between 303 and 316 rounds, 
at an average of 309.0 ± 2.5 rounds. With ksdyn = 1.0 and Sadd = 10 it took with 
Manhattan metric 176.9 ± 1.8, Chebyshev metric 235.2 ±2.2, and VI metric 173.4 ± 1.7. 



6.3. A Crowd Moving Along a Straight Corridor 

As the dynamic distance potential field method can lead to an increase of the walking 
speed, one has to show that the reduction of the evacuation time as exhibited in table 
|3]is a consequence of the agents utilizing space at the corner more efficiently and not a 
consequence of them just walking faster. Therefore simulations have been done, where he 
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Figure 6. Situation after 90 seconds with the dynamic distance potential field method. 
The ordering of figures with regard to different methods and parameters is the same 
as in table [3l 



agents did not have to move around the corner but walk on straight forward. Anjd;hing 
else was left as shown in figure |H Without dynamic distance potential field the process 
took between 113 and 118 rounds, at an average of 113.5 ± 1.6 rounds. The results for 
the various metrics and different parameters is shown in HI 

The results from H] show two remarkable things: first off, the dynamic distance 
potential field does not speed up the process. Second, the combination of Manhattan 
metric with large ksdyn can - at least in this case with motion along the axis of 
discretization - lead to considerable delays. The reason is not that each single agent 
is slowed down directly, but that small density variations build up as the agents 
simultaneously move to lanes where in front of them relatively few other agents are 
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ksdyn I Sadd 



Manhattan 



Chebyshev 



VI 



1.0 / 2 
3.0 / 2 
1.0 / 4 
3.0 / 4 
1.0 / 10 
10.0 / 10 



116.6 ± 1.4 (114, 121) 

133.0 ± 1.6 (130, 139) 

118.1 ± 1.5 (116, 122) 
137.3 ± 1.5 (135, 144) 
119.3 ± 1.5 (117, 125) 
154.1 ± 1.9 (151, 163) 



113.8 ± 1.2 (112, 117) 
114.8 ± 1.7 (112, 122) 
114.1 ± 1.5 (112, 120) 

115.3 ± 1.6 (113, 121) 

114.4 ± 1.4 (112, 119) 
115.8 ± 1.5 (113, 120) 



114.1 ± 1.5 (111, 118) 

114.6 ± 1.3 (112, 118) 

113.8 ± 1.2 (111, 117) 

114.9 ± 1.4 (112, 118) 

113.7 ± 1.3 (111, 117) 
116.6 ± 1.5 (114, 123) 



Table 4. Number of rounds until the last agent has walked through the corridor and 
left the simulation 



located. This position is evaluated as being particularly bad and the agents move to 
the side in subsequent rounds. But it's not the sidewards motion that leads to delays 
but the locally high density right in front of them toward the destination. Figure [7] 
illustrates this effect. As this under any circumstances can only be an unwanted effect 
the use of the Manhattan metric demands greater care in setting ksdyn and after the 
simulation is carried out a check, if or if not this artificial lane formation appeared at 
some time. It is therefore better to not use the Manhattan metric at all or to replace 
ksdyn by another parameter a from which ksdyn calculates for example as 

ksdyn = (10) 

In this case ksdyn can only be smaller than ks- 

The Chebyshev metric on the contrary only has a minor influence on both the 
shape of a crowd moving around a corner and the evacuation time. What remains is 
the method of Variant 1. Nevertheless for the rest of the contribution consider all three 
metrics will be considered to check and consolidate the impression that there are good 
reasons to use variant 1 (VI). 



Figure 7. The situation in three subsequent rounds: 89 (red), 90 (green), and 91 
(blue) when Manhattan metric is used. It can be seen that there are specific lanes in 
each round and that the lanes of immediately consequent rounds are almost entirely 
separated (there is no yellow or cyan), but that there is one lane in round 91 at a 
position where there was one in round 89 (magenta). 



6.4- A Crowd Doing a U-Turn Around a Wall 

Next the agents from the same starting positions had to do a u-turn around a 40 cm wide 
wall. Compared to figure H] the unaccessible (white) area above the starting corridor 
were accessible and the destination was on the other side of the wall "above" the flow 
restricting columns. Without dynamic distance potential field the process took between 
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461 and 476 rounds, at an average of 468.75 ±3.1 rounds. The results for the various 
metrics and different parameters are shown in [51 



ksdyn I Sgdd 



Manhattan Chebyshev VI 



1.0 / 2 335.7 ± 2.1 (330, 341) 385.9 ± 2.1 (382, 392) 375.3 ± 2.7 (370, 384) 

3.0 / 2 273.5 ± 2.1 (268, 279) 331.7 ± 2.2 (325, 338) 313.2 ± 2.2 (309, 322) 

1.0 / 4 260.5 ± 1.8 (255, 265) 334.9 ± 2.3 (329, 342) 299.5 ± 2.1 (294, 304) 

3.0 / 4 234.9 ± 2.0 (231, 240) 302.1 ± 2.1 (297, 308) 249.1 ± 1.8 (245, 255) 

1.0 / 10 225.0 ± 2.1 (219, 231) 315.4 ± 2.4 (310, 320) 249.8 ± 1.8 (246, 256) 

10.0 / 10 233.0 ± 2.4 (229, 239) 298.5 ± 3.2 (291, 307) 235.1 ± 2.4 (228, 243) 

Table 5. Number of rounds until the last agent has done the u-turn and left the 
simulation 

These results reveal a reduction of the evacuation time by 50% for some of the 
parameters compared to the case without dynamic distance potential field. Again the 
Chebyshev metric has the smallest and the Manhattan metric the largest effect. The 
strength of artificial effects strongly depends on ksdyn and almost not on Sadd- 

6. 5. The Fundamental Diagram 



In 28| it was shown that the F.A.S.T. model can reproduce Weidmann's fundamental 



diagram [M] quite well with the rather simple set of parameters Vmax = 4, fc^ = 1.2, and 
kother = 0.0. Therefore the fundamental diagram was calculated with these parameters 
and the djTiamic distance potential field added. And it was calculated using the same 
geometry: a 4 meters (10 cells) wide ring shaped corridor with an outer radius of 500 
cells. There were four destinations at 0, 90, 180, and 270 degree making the agents 
walk around. The reason for the selection of this geometry once was that it would show 
remaining anisotropics and that in principle one could implement this in an experiment, 
which is not possible in a straight corridor with periodic boundary conditions. The 
results are shown in figure [H 

There are two possible reasons for the reduction of the flow at the intermediate 
densities: a) the deviation of the dynamic metric from the Euclidean metric and b) 



effects of percolation [45[ . To decide about this an additional fundamental diagram was 
calculated with four additional intermediate distances at 45, 135, 225, and 315 degree 
and VI metric. As VI produces errors only at corners, with these additional intermediate 
destinations, the average error of the dynamic distance potential field was reduced. In 
effect the flow actually increased, but the maximum effect was an increase by 0.02 agents 
per round and cell at a density of 0.3 agents per cell. So, by far the largest share of flow 
reduction must be due to the percolation of unoccupied cells vanishing with increasing 
density. This is supported by the fact that both for Manhattan and Chebyshev metric 
the largest deviations from the case without dynamic distance potential field appear at 
densities about 0.06 agents per cell smaller than the percolation threshold of agents (on 
an infinite lattice) with regard to the metrics' neighborhood. Analytical investigations 
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0,5 




Density [agents / cell] 

Figure 8. Theoretical fundamental diagram (fd) by Weidmann (see equation ((TT|)), 
fd from the F.A.S.T. model without dynamic distance potential field, and with the 
dynamic distance potential fields using Manhattan, Chebyshev, and VI metric. All 
F.A.S.T. simulations have been carried out with ks = 1.2, those with dynamic distance 
potential field with ksdyn = 1-0 and Sadd = 10. All other k were set to zero. As one 
round equals one second and one cell has a size of 0.4 X 0.4 sqm, the density is shown 
in fractions of 6.25 agents per sqm, and the specific flux in fractions of 2.5 agents per 
meter and second. 



of this phenomenon are difficult not only due to the geometry, but as one cannot assume 
in advance that the agents' positions are uncorrelated. 

Annotation: Weidmann's theoretical fundamental diagram follows the formula 

with flow j, density p, free speed Vfree = 1-34 m/s, maximum density pmax = 5.4 1/sqm, 
and the gauge constant 7 = 1.913 1/sqm. 

6. 6. A Part of a Stadium 

To go beyond the elementary geometries investigated so far a scenario with 10.000 agents 
moving through a simplified stadium-like geometry was simulated. Figure [9] shows the 
geometry and initial positions of the agents. The average time of 100 runs until the last 
agent had left the scenario without the dynamic distance potential field method was 
2956.8 ± 8.1 rounds (seconds). The average of individual egress times was 1579.4 ± 3.5 
rounds. The results for the process using the various metrics and different parameters 
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Figure 9. Geometry and agents' starting positions for the simplified stadium test 
simulation. On the magenta and dark green areas two times 5,000 agents started their 
egress and had to walk to the light green area on the upper right. 



ksdyn / Sadd 


Manhattan 


Chebyshev 


VI 


1.0 / 2 


2096.2 


2228.2 


2040.4 


3.0 / 2 


1785.9 


1900.1 


1697.4 


10.0 / 2 


1739.6 


1854.3 


1753.8 


1.0 / 4 


1673.9 


1905.9 


1620.6 


3.0 / 4 


1568.1 


1823.8 


1424.4 


1.0 / 10 


1503.5 


1861.5 


1417.5 


10.0 / 10 


1646.5 


1959.0 


1962.7 



Table 6. Average of 100 number of rounds until the last agent has left the simplified 
stadium simulation. 

The first fact to reahze when looking at tables [6] and [7] is that method VI leads 
to the smallest total and individual evacuation times for ksdyn = 1-0 and ksdyn = 3.0. 
With VI the evacuation times are larger for ksdyn = 10.0 than for the next smaller ksdyn 
with identical Sadd- In total the effects show most pronounced - and also "best" - for 
Sadd = 10 and unrealistic artifacts are smallest for ksdyn = l-O. Therefore in figures 
[To] to [19] the case without dynamic distance potential field is compared at t = 400 and 
t = 800 to the three variants using this choice of parameters. It's interesting to take a 
look at 

• The shape of the crowd waiting on the initial area. 
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f^Sdyn 1 ^add 


IVTpi Ti n f f PTi 

IV J.CXX111CX u v CXiLl 


(! n phvc; h pa/" 


VI 

V X 


1 n / 9 


1 147 






3.0 / 2 


1003.2 


1084.8 


980.5 


10.0 / 2 


992.0 


1067.4 


1045.0 


1.0 / 4 


950.6 


1094.8 


936.4 


3.0 / 4 


916.6 


1055.2 


856.2 


1.0 / 10 


886.2 


1077.6 


842.8 


10.0 / 10 


957.3 


1109.1 


1119.7 



Table 7. Average of 100 individual egress times in the simplified stadium simulation. 

• The vertical position of the stream walking from left to right. 

• The jam at the corner, where the main moving direction changes. 

• The distribution of agents in the final corridor. 

• The utilization of the gates at the end. 

For all of these issues things are either looking best using VI or there is no difference 
between the three metric variants. In any case it's better to make use of a dynamic 
distance potential field. 

Apart from the simulations for which the results are shown in tables [6] and [7] the 
following variants have been simulated: 

• All three metrics with ksdyn = 1-0 and Sadd = 100. The results did not in any aspect 
(evacuation time, individual egress time, shape of crowd, etc.) deviate significantly 
or recognizably from the ones with Sadd = 10. This confirms that with increasing 
Sadd there is a saturation of the influence of this parameter. 

• VI again with ksdyn = 1-0 and Sadd = 100, but now the field values for Sdyn and S^^y^ 
were not rounded to integer values after taking the square root in the calculation 
process. This just as well did not create recognizable macroscopic effects. 

• VI again with ksdyn = 1-0 and Sadd = 100, but now the ASdyn values for the 
current position of a agent were ignored (artificially set to zero). This created a 
very unrealistic behavior, long gent had other agents in some not too 
large distance ahead of him he preferred his current position a lot over any other 
and remained standing still. 

• Without dynamic distance potential field, but with kw > and kj > 0, so the 
static repulsion from walls and inertia were switched on. With that the evacuation 
time could be reduced. With k^ > but kj = there was only a minor effect of a 
reduction of 30 rounds at kw = 0.5. From all parameter combinations investigated 
the largest effect was seen with k^ = and fc/ = 2.0 as the average evacuation 
time reduced to 2026.5 rounds. However, this cannot be attributed to the jamming 
crowd at the corner making following agents detour, but to a wider spread of the 
stream of agents coming from the upper source due to inertia. This spread would 
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have taken place just as well, if there had not been a jam at all at the corner. In 
this case the increased inertia would have led to an increase of evacuation time. 
Thus, the behavior induced by the parameters was independent of the situation 
and therefore only accidentally led to a reduction of evacuation time. 




Figure 10. The situation after 400 rounds without dynamic distance potential. 




Figure 11. The situation after 400 rounds with Manhattan metric. 
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Figure 12. The situation after 400 rounds with Chebyshev metric. 




Figure 13. The situation after 400 rounds with VI metric. 



6.7. Computation Times 

The stadium was simulated for all three metrics with Sadd = 1? taking the effect from 
the method. The computation time is compared to the average simulation time of 6:00 
with the calculation for the method switched off. The results are shown in table [8] 

When looking at the computation times during the calculation of the fundamental 
diagram, one can see the difference in computation time at different densities, 
respectively different numbers of agents in identical geometries. The computation times 
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Figure 14. The situation after 800 rounds without dynamic distance potential. 




Figure 15. The situation after 800 rounds with Manhattan metric. 





Manhattan 


Chebyshev 


VI 


with dyn. potential 


6:30 


6:48 


7:18 


relative increase 


8.3% 


13.3% 


21.6% 



Table 8. Average computation time ([minisec]) with the dynamic potential field 
method with Sadd = 1 and the relative increase compared to simulations without 
dynamic distance potential field method. 
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Figure 16. The situation after 800 rounds with Chebyshev metric. 




Figure 17. The situation after 800 rounds with VI metric. 



in figure [20] were obtained with Sadd = 1, so the behavior of the agents did not change 
compared to the case without dynamic distance potential field. It can be seen that the 
absolute difference is about constant, while the relative one decreases with increasing 
number of agents. 

The relative increase in computation time depends on the scenario, as the time 
for the calculation of the dynamic distance potential field only depends on the size of 
the area, but almost not on the number of agents in the simulation. This is a very 
positive fact compared to a naive ansatz of each agent having to observe any other to 
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Figure 19. The situation after 800 rounds with Manhattan (red), Chebyshev (green), 
and VI (blue) metric. 




100 

-I , , , , , 

5000 10000 15000 20000 25000 30000 

^Agents 

Figure 20. Average computation times (in seconds) for various densities in the 
fundamental diagram. 1100 rounds were simulated. 

meet the goals as presented here. With such a naive ansatz the computation time would 
depend on the square of the number of agents. With the dynamic distance potential 
field the relative increase of computation time decreases with increasing global density, 
as typically the computation time depends at least linearly on the number of agents. 
As the method is most valuable to make agents evade large aggregations of other agents 
(jams) early, a user can switch off the method in simulations, where the local density is 
small everywhere and all the time. So, in the case when the method would imply the 
largest relative computation time increase, it is of no use anyway and can be skipped in 
the simulation algorithm. 

With regard to usage in other models, the relative increase in computation time is 
small, if the model needs much computation time. The computation time of the dynamic 
distance potential field is model independent. If a model has many time steps to simulate 
one second it might be possible to recalculate the dynamic distance potential field only 
once a simulation second, as within a second the situation never changes dramatically 
in agent dynamics. So, if there is a model with the same computational effort per time 
step as the F.A.S.T. model, but which is doing 20 time steps per simulation second, the 
increase in computation time might only be 1%. 
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7. Conclusions and Outlook 

The method proposed in this work proofed to do what was intended: the movement of 
the agents looked more natural in many cases. The agents showed more jam avoiding 
behavior, accepting longer paths around corners in exchange for higher walking speeds 
and thus the group of agents increased the total flow. The effect was strongest for the 
Manhattan metric. The smallest effect was achieved using the Chebyshev metric. Partly 
this may be due to the alignment of the geometry of most examples of this article to the 
axis of discretization. Without empirical check, it seems that the shape of the crowd 
is most natural using the VI metric. As it is a combination of the two other methods 
effects on S^yn can reach far, but decrease with distance. It seems that setting Sadd to 
large values close to saturation (such that further increase does not change anything 
anymore) pronounces the positive effects most without bringing in negative drawbacks. 
It is possible, however, to choose too large values for ksdyn- Then the crowds will form 
unnatural shapes. In situations where the introduction of the method would not have 
been necessary the method only introduces small changes - at least if ksdyn is reasonably 
small. 

A metric investigated in [32I, but not considered here, is V2, where the Moore 
neighborhood is used but a is added for in the flood fill over corners. As a real- 
valued VI metric seems to bring no changes, the integer-valuedness of the metric seems 
to be not necessary and V2 can also be a candidate for investigation. This would more 
or less only be formally, as VI works well enough. It might be, though, that V2 works 
just as well and is quicker. For such an investigation there would be two variants. In 
one Sadd would be added in steps over edges and corners alike, if the corresponding cell 
is occupied. In the other method one would multiply 1 or \/2 with a parameter Smuit- 

In certain situations the influence of large jams progresses over long ranges, 
sometimes infinitely far. This was least pronounced for Chebyshev metric. While 
the movement of the agents at first seems to be unrealistic using Chebyshev metric, 
as the limited range of influence makes the agents try to evade the jam only much 
shorter in front of the jam, this could be used to simulate the influence of obstacles, 
smoke or darkness, where visibility range is limited as well. However, the trajectories 
generated by the influence of Chebyshev dynamic distance potential fields only create 
some associations of the behavior of real pedestrians in such situations. It might well 
be that actually there are large differences. 



In [32| it was shown that the error of VI depends on the number of corners the 
potential has to flow around. Scenarios with a large number of such corners were not 
investigated in the present work. So it might be that for such geometries the VI metric 
cannot be used as it might produce unrealistic results. This might be hard to prove, 
though, as the complexity of such a scenario a) makes it hard to develop an idea of how 
real pedestrians might walk and b) if this is possible, it might still be difficult to clearly 
detect deviations from this behavior in the simulation. 

Whether one uses the method proposed in this contribution or not: a result of 
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repeated simulation of pedestrian dynamics is the field Ti{x,t) of travel times from 
some position to the destination. Just like the correction of the static distance potential 
field by the dynamic one, Ti (x, t) is or at least should be a better estimation of the 
real travel times. Therefore one could use Ti{x,t) similar to Sdyn{x,t) in a second set 
of simulations to calculate a T2{x,t). The hope is that simulating iteratively converges: 
|Tj+i(x, t) — Ti{x,t)\ < e V z > j for a fixed j. This would be a generalized dynamic 
assignment [46] for fully two dimensional application, similar as in For such an 

iterative strategy one would have to define an extrapolation method for the travel times 
of those locations that were not accessed by agents in the previous iteration. The 
importance and the requirements for the quality of this method are the higher the 
more the initial simulations deviate from the final solution. In this respect the method 
proposed in this work can be assumed to be very helpful. 
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